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Abstract 

We extend our model for the pion, which we used previously to calculate its diagonal structure function, to the off-forward 
case. The imaginary part of the off-forward 7*7r — > 7*7r scattering amplitude is evaluated in the chiral limit (m^ — 0) and 
related to the twist-two and twist-three generalised parton distributions H, H 3 , H 3 . Non-perturbative effects, linked to the 
size of the pion and still preserving gauge invariance, are included. Remarkable new relations between H, H 3 and H 3 are 
obtained and discussed. 
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1 Introduction 

Structure functions, which can be extracted from deep- 
inelastic experiments, are useful tools to understand the 
structure of hadrons. Even if their Q 2 evolution is con- 
sistent with perturbative QCD, they result mainly from 
non-perturbative effects that are still not calculable in 
the framework of QCD. This has led to phenomenolog- 
ical quark models embodying various non-perturbative 
aspects of QCD. These models can be used to depict the 
behaviour of the structure functions and to understand 
the connection between data and non-perturbative as- 
pects of hadrons. There has been extensive work on di- 
agonal distributions along these lines (see Refs [l]-[2] for 
the pion case). These distributions can be used as the 
initial condition for a DGLAP evolution, which is neces- 
sary before a comparison with data [3] . Such models can 
be applied to the off-diagonal case, for which generalised 
structure functions [4] can be linked [5] to generalised 
parton distributions (GPD's). 

One of the setbacks of phenomenological quark models 
suited to the description of the low-energy features of 
hadrons is that the underlying quark structure is ob- 
scured by the necessary introduction of regularisation 
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procedures which result in non-negligible differences in 
the structure functions. 

To avoid these complications, we investigated the diag- 
onal structure functions in the case of a simple model 
for the pion [1], where the pion-quark-antiquark pseu- 
doscalar coupling {rfg-nqq) yields the correct symmetry, 
while the non-perturbative aspects come from a momen- 
tum cut-off mimicking the size of the pion, but still pre- 
serving gauge invariance. This freed us from the ques- 
tion of what would be the detailed inner structure of the 
meson. In that calculation, owing to the introduction 
of such a cut-off, crossed diagrams for the pion-photon 
scattering appear as higher twists, leading twist struc- 
ture functions can be identified, and a reduction of the 
momentum fraction carried by the quarks is observed. Of 
course, as the cut-off is relaxed to let the quarks behave 
freely, the momentum sum rule (2x) = 1 is recovered at 
infinite Q 2 . Having that tool at hand, we now turn to 
the investigation of the properties of off-diagonal par- 
ton distributions, which are likely to shed some light on 
parton correlations and which have therefore attracted 
much interest in recent years [6,7,8,9,10,11]. 

In the following, we calculate the imaginary part of the 
off- forward photon-pion scattering amplitude, and of the 
structure functions F\ , . . . , F$, related to the five inde- 
pendent tensor structures in the scattering amplitude, 
and we discuss their behaviour. We relate them to vector 
and axial vector form factors and to the twist-two and 
twist-three generalised parton distributions (GPD's) H, 
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H 3 and H 3 . We shall show that, within our model and in 
the high-Q 2 limit, the non-diagonal structure functions 
F 3 and F4 are related to Fi, while F 5 happens to be a 
higher twist. These results lead to new relations for the 
GPD's in the neutral pion case. 



have included a factor 1 /m 2 in the definition of F$ in or- 
der to avoid divergences when the chiral limit is taken. 
Note that Bose symmetry requires f\, F 2 , F 4 , F 5 to be 
even and F3 to be odd in £. 



2 General tensorial structure of the 7*7r — > 7*77 
amplitude 

2. 1 External Kinematics 

Let pi and P2 be the momenta of the ingoing and out- 
going pions, qi and q2 those of the corresponding pho- 
tons (see Fig. 1). Defining p = (pi + P2) , q = an d 
A = P2 — Pi , one can then write the scattering ampli- 
tude as a function of the Lorentz invariants t — A 2 , 



Q 2 



-q 2 , x 



Q 

2p-q 



and £ 



In the elastic limit, characterised by (g+^) 2 = (g— ^) 2 , 
one has A • q = and thus £ = 0, while the diagonal 
limit (A = 0) is obtained for £ = and t = 0. We 
further recover the Bjorkcn variable x — xb 1 where = 
il_ 

2pi-<?i ' 

For Virtual Compton Scattering (VCS), for which the 
outgoing photon is on-shell, £ is related to x through £ = 

— x ^1 — • Hence in the Deeply Virtual Compton 
Scattering (DVCS) limit, t < Q 2 and £ = -x. 



2.2 The structure functions Fi 's 
The hadronic tensor is defined through 

T^(q,p, A) = i I dVc-s (p 2 |Tj M (r/2)j I ,(-r/2)|p 1 ) .(1) 



There exist five independent kinematical structures in 
Eq. (1) that parametrise the photon-pion amplitude. 
Defining the projector = — 9 ^ 9 ^ and making 
use of these five structures, we can rewrite as follows: 



T^(q,p,A) = -V^g aT V Tl ,F 1 



r^P a p T Vrv 
p ■ q 



H ~ ^3 



+ ■ 



7V(p CT (A- - 2£p T ) - (A- - 2fr> T )7> TI , 



2p-g 

+^ CT (A ff - 2CK)(A T - 2fr T )V TV F b . 



(2) 



Current conservation is ensured by means of the projec- 
tor V,j, v . Our notation slightly differs from Ref. [5]: we 



3 The model 

3.1 General description 

We use the pion model introduced in our previous 
work [1], in which the qqn vertex is represented by 
the simplest pseudoscalar coupling. The Lagrangian 
includes massive pion and massive quark fields interact- 
ing through the pseudoscalar vertex, with an effective 
pion-quark coupling constant. 

Considering an isospin triplet pion field 7? = (tt + , 7r° , 7r~ ) 
interacting with quark fields tp the Lagrangian density 
reads 



where f is the isospin vector operator. 



(3) 



Of course, if our pseudoscalar field is to represent real pi- 
ons, we have to impose that the corresponding hadrons 
have a finite size. That we shall do through the use of 
a cut-off, as detailed below, the choice of which sets a 
constraint on the value of the quark-pion coupling con- 
stant [1]. 

We shall limit ourselves in this paper to the calculation 
of the imaginary part of the scattering amplitude, which 
allows a direct comparison with our previous work and 
which is sufficient to determine the GPD's of neutral 
pions [5]. 

At the leading order in the loop expansion, four diagrams 
contribute. They are displayed in Fig. 1. Following the 
kinematics defined in section 2.1 and applying Feynman 
rules, it is straightforward to write down the analytical 
expression for the scattering amplitude. For a given set 
/i, v of the photon indices and with well-known conven- 
tions 4 , the contribution of the first diagram (a) shown 
in Fig. 1 to the scattering amplitude reads 



The isospin/charge factor (e„ + e 2 d ) corresponds to the 



following choice of the isospin matrix: it 
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V2 0, 
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Fig. 1. The simplest diagrams contributing to the imaginary part of the amplitude for the scattering 7*71 — > Upper (lower) 

diagrams are referred to as box (crossed) diagrams. Dashed lines represent the discontinuity of the amplitudes, i.e. their imaginary 
parts. 
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d 4 fcTr( 7 5 ((|(;-^) + m g )7 5 
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r).(4) 



According to Cutkosky rules, the imaginary part of 
the amplitude is obtained by putting the intermediate 
quark lines on shell. This is realised by the introduc- 
tion of the two delta functions, 5 ((fc + q) 2 — m 2 ) and 



Working out the delta functions, we obtain that: 



Expressions from the other three diagrams of Fig. 1 - the 
one with reverse loop-momentum and the two crossed 
diagrams - are similar and are not written down. Results 
below all pertain to the chiral limit m-^ = 0. 



5((k + q) 2 -m 2 q )8((k-p) 2 -ml) = 

— —r S (cos 9 - cos 6 ) ■=- r 6 (fed - k' ) 
2k p \q\ 2y/s 



(7) 



3.2 The implementation of the cut-off 

A simple way to impose that the pion has a finite size is to 
require that the square of the relative four-momentum of 
the quarks inside the pion is limited to a maximum value 
A 2 . Before writing this explicitly, let us give the details of 
the internal kinematics, i.e. the one involving the loop- 
momentum k. Let (j) and be the spherical angles of 
k with respect to the z-axis taken as the direction of 
the incoming photon. Defining k 2 p = \k\ 2 and r = k 2 = 
k 2 



k 2 and using spherical coordinates, we write the 



element of integration as: 
d A k = dkodk p k 2 d(cos6)d4>: 
or with the help of the variable r: 
k 

d A k — dko-^drd(cos 



(5) 



(0) 



with cos#o 



2k qo-Q — m +t 



and fcn 



Fi- 



2k p \q\ '"^ n 2v ^ 

nally, the element of integration over the internal mo- 
mentum, considering only the imaginary part of the am- 
plitude, reads: 



d k = drdtfi 



kQ—k' Q ,cos 9— cos f?o 



(8) 



with \q\ = ^L ^s'Q'+a-zs) 2 ^ . The boundary values 
of the integration domain on r are obtained by solving 
cos $0 = ±1- 



Now we may look at the effect of the finite size of the 
pion on the integration procedure upon k. The relative 
four-momentum squared of the quarks inside the pion is 
given by 
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Of=(2k-p±^j 

= 2t + 2m\ - ml + l - ± 2k • A, (9) 

for pion-quark vertices like the ones in diagram 
Fig. l.(a), and by 

0± - ^2k-p + 2q±^j 

= -2T + 6m 2 q -ml+ t --^±2(k-A + ^l), (10) 

y 2 x x 

for vertices as in diagram Fig. l.(b). Note that k ■ A is a 
known function of the external variables as well as of 9 
and r. Generalising the procedure of [1], we require ei- 
ther | Of | < A 2 or \Of\ < A 2 for all dia grams. Gauge 
invariance is preserved by the cut-off, since \Of\ depend 
only upon the external variables of the 7*7r — > qq pro- 
cess. 

As the Oj's and r are always negative, we require one of 
the two following conditions: 



J F 1 {x)dx = A (12) 



As Fi/g 2 is a priori a function of Q 2 , the sum rule im- 
poses that g should be a function of Q 2 . But at high 
enough Q 2 , where the details of the non-perturbative 
interaction are less and less relevant, Fi/g 2 reaches its 
asymptotic shape when the cut-off procedure is ap- 
plied, and we obtain a constant value for g. In Ref. [1], 
this asymptotic regime was reached for Q 2 as small as 
2 GeV 2 . In the following, we shall make use of these 
previously obtained values, which are functions of the 
cut-off A. 

However in the DVCS case, an ambiguity may arise as 
one of the vertices has an external kinematics similar to a 
vanishing Q 2 DIS. This ambiguity is lifted if one notices 
that the pertinent quantities are not q\ and q\ separately 
but the factorisation scale, which may be taken as the 
square of their mean, Q 2 . Thus in DVCS, although q 2 
vanishes, Q 2 does not and we shall consider that g is 
constant. 

4 Results 

4 ■ 1 General features 



-A 2 mi 2 t ., .. 
r> — + -f-m 2 q -- + \k-A\, 

For t small, |Oi| and | O2 1 cannot be small simultane- 
ously. The crossed diagrams have their main contribu- 
tion for 0\ ~ O2, and are thus suppressed by a power 
A 2 /Q 2 when the cut-off is imposed. The box diagrams 
have a leading contribution for or | O2 1 small, and 
are not power suppressed by the cut-off. 

It may be worth pointing out that the vertical propaga- 
tors are more off-shell in DVCS than in DIS, hence one 
would expect DVCS to be better described by perturba- 
tion theory than DIS. 

3.3 The coupling constant 

In the diagonal case, we have determined the coupling 
constant g — g^ qq by imposing that there are only two 
constituents in the pion. This sum rule constraints F\ as 
follows [1]: 



From the imaginary part of the total amplitude, the five 
structure functions Fi can be obtained by a projection on 
the corresponding tensors. From now on, Fi will stand for 
the imaginary part of these structure functions. In order 
to display their general features, we plot them in Fig. 2 
first as functions of x and £ for parameter values m q — 
0.3 GeV and A = 0.75 GeV, to ease the comparison 
with [1], and for Q 2 = 10 GeV 2 and t = -0.1 GeV 2 . 
Let us notice that for any fixed value of £ not close to 
±1, we recover for F\ and F 2 the same behaviour as in 
the diagonal case. We checked indeed that the diagonal 
limit is recovered for £ = and t — 0. Let us notice also 
that the structure functions F 3 , F4, F 5 depend little on 
£ except when this variable is close to ±1. 

Let us turn now to DVCS. Fig. 3 displays the behaviour 
of Fi for various values of t with and without cut-off. In 
the presence of size effects, the value of F\ gets signif- 
icantly reduced, especially for small x, as \t\ increases, 
whereas that effect is much less noticeable without cut- 
off. 

In the elastic case (see Fig. 4), the same suppression 
at small x is observed, especially when the cut-off is 
applied. In Fig. 5, we display the average value of 2x with 
respect to the F\ distribution. The value of (2x) increases 
when \t\ increases. The momentum fraction carried by 
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Fig. 3. F\ as a function of x for various values oft in the DVCS case with (a) and without (b) cut-off (A = 0.75 GeV), at 
Q 2 = 10 GeV 2 . 
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Fig. 4- Evolution of Fi (elastic case, £ = 0) for decreasing values oft with (a) and without (b) cut-off, for Q 2 = 10 GeV 2 
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Fig. 5. Mean value of2x for the Fi distribution in the elastic case (£ = 0) as a function of Q 2 and for various values oft. 
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Fig. 6. F\ as a function of x fort — —0.1 GeV 2 and for various values of Q 2 in the DVCS case with (a) and without (b) cut-off. 



the quarks and probed by the process thus increases with 
the momentum transfer. 

The effects of the variation of Q 2 are displayed in Fig. 6. 
As in the diagonal case [1], we can conclude that the 
details of the non-perturbative effects cease to matter for 
Q 2 greater than 2 GeV 2 , that is significantly larger than 
A 2 . On the other hand, when the cut-off is not applied, 
we sec (Fig. 6 (b)) that F\ evolves so slowly with Q 2 that 
the asymptotic state is not visible. 



4-2 High-Q 2 limit: new relations 

Having determined the 5 functions F^s in the context of 
our model, we shall now consider their behaviour at high 
Q 2 . Expanding the ratios of ,|f , Jf , we obtain the 

following asymptotic behaviour: 



F 2 = 2xFi + 0(l/Q 2 ) 

^3 



f, = ^Kf 1 + o(i/q 2 ), 



F 4 = ^F 1+ 0(1/Q 2 ) 
F 5 =0(l/Q 2 ). 



(13) 
(14) 

(15) 
(16) 



The fact that at leading order there are only three in- 
dependent structure functions has been known for some 
time [5,10]. However, we show here that they can all be 
obtained from F\. The first relation is similar (at lead- 
ing order in 1/Q 2 and with the replacement of x by xb) 
to the Callan-Gross relation between the diagonal struc- 
ture functions F\ and F 2 , valid for spin one-half con- 
stituents in general. Except for F$, which is small at 
large Q 2 , these relations show that F 2 , F3 and F4 are 
simply related to F\ at leading order. These relations 
clearly display and therefore confirm the symmetries of 



these functions. Combining Eqs. (14) and (15), we have, 
at leading order, 



^3 = ^4, 



(17) 



which confirms that F 3 is an odd function of £, while F4 
is even. 

The simple relations between the Fj's (at leading order) 
constitute a remarkable result of our model. Further- 
more, we checked that the term 0(1/ Q 2 ) in Eq. (13) is 
numerically quite small, even for moderate Q 2 . One may 
wonder whether these results are typical of our model or 
more general. 



5 Linking the F 4 's to H, H 3 , and H 3 

Having at hand the five functions Fi 's that parametrise 
the amplitude for 7*7r — ► 7*7r, we would like to link 
them to the off-forward parton distribution functions or 
to the generalised parton distributions. For this purpose, 
we make use of a tensorial expression coming from the 
twist-three analysis of the process, which singles out the 
twist- two Ti and the twist-three W 3 ,7l! 3 form factors. 
Following Ref . [5] , we write 5 : 



q-Vi 
~2p-q 

'2p-q 



V 

p ■ q 



(18) 



5 Please note that Ref. [5] uses V VI1 instead of V^v as pro- 
jector. 



7 



"tq 



: 0.1 

■- 0.3 
■- 0.7 
■- 0.9 



0.5 



-0.5 



0.2 0.4 0.6 0.8 



1 


e = o.i 




= 0.3 




= 0.7 




= 0.9 


1 


i i i 



0.2 0.4 0.6 0.8 




Fig. 7. A comparison of the values of Hz and Hz obtained in our model (left) with those calculated in the Wandzura-Wilczek 
approximation, using the value of H from our model (right), for Q 2 = 10 GeV 2 . 



where the V^'s and A\ read 
Vi P = 2 Pp H + (A p - 2& p )H 3 + twist 4, 
Aip 



ItpApq ^3 



p-q 



V, p = xV lp -^ q -V 1+ l - e -^-Al 
2 p ■ q 4 p ■ q 



twist 4. 



(19) 
(20) 

(21) 



In Ref. [5], gauge invariance of Eq. (18) beyond the twist- 
three accuracy was in fact restored by hand, contrarily 
to the present calculation for which the amplitude is 
explicitly gauge invariant. 

To relate the F^s to the Ws, we project the ampli- 
tude (18) onto the five projectors contained in Eq. (2) 
and identify the results with the i^'s. Note that, in the 
neutral pion case, the imaginary part of the form factors 
H directly gives the GPD's H, H 3 and H 3 up to a factor 
2ir. As we have kept the off-shellnesses of the photons 
arbitrary, we in fact can relate the imaginary parts of Fi 
to the GPD's for arbitrary x and £: 



1 

27 
1 

27 

2lT X 



Fx=H, 



F 2 = 2xH + 0{-^), 



i 

Q 2 



(22) 
(23) 
(24) 



-F 4 = -— 2 (H^x + H 3 x 2 Hx) + O(-), (25) 



(26) 



Replacing the i^'s by the expressions (13-16), we can 
write 



H 



X (e - 1) 



(27) 



and 



"^iw^^ ^-' 1 ^ ^- (28) 



As Fi to F4 can be written in term of only one of them, 
e.g. Fx, it is not surprising that H 3 and H 3 are simply 
related to H. Note that polynomiality of the Mcllin mo- 
ments of H, H 3 and H 3 , together with Eqs. (27) and 
(28), imply that H must be a polynomial Pjj multiply- 
ing £ 2 — 1. The fact that, as can be seen from Fig. 7, H3 
is almost independent of £ shows that Pr is very close 
to a constant. 

To convince ourselves that relations (27) and (28) are 
new, we have compared them to the Wandzura-Wilczek 
approximation [12], given in the pion case in [5,13]. First 
of all, it is well-known that these relations are discon- 
tinuous at £ = ±x, which is not the case for (27) and 
(28). Furthermore, we show in Fig. 7 the results of the 
Wandzura-Wilczek approximation compared with our 
results. We see that the two are numerically very differ- 
ent. Hence the relations (27) and (28), derived in an ex- 
plicitly gauge- invariant model, do not come from " kinc- 
matical" twist corrections, but emerge from the dynam- 
ics of the spectator quark propagator and from finite-size 
effects. 
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6 Discussion and conclusion 

We have extended our previous model for the pion to 
investigate the off-diagonal structure functions for this 
particular case. The introduction of a cut-off allows the 
crossed diagrams to behave as higher-twists and to relate 
the imaginary part of the forward amplitude with quark 
GPD's. 

We used the formalism of Ref. [5] in order to decompose 
the amplitude along the relevant Lorentz tensors, to de- 
fine five structure functions Fi , and to relate the latter to 
the GPD's H, H 3 and H 3 introduced in the twist anal- 
ysis. We have found that our results in the forward case 
are qualitatively preserved when departing from the for- 
ward limit. 

Our investigation yields new results. In particular, we 
singled out new relations, which link the F^s in a simple 
manner at leading order in 1/Q 2 . More intriguing, we 
found that the twist-three structure functions are simply 
related to H by relations that differ from the Wandzura- 
Wilczck approximation. 

Although these relations are derived in the context of our 
simple model, it is possible that they can be extended 
to a more general case. 
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